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Certain dissipative Ginzburg-Landau models predict existence of planar 
interfaces moving with constant velocity. In most cases the interface solutions 
are hard to obtain because pertinent evolution equations are nonlinear. We 
present a systematic perturbative expansion which allows us to compute 
effects of small terms added to the free energy functional of a soluble model. 
As an example, we take the exactly soluble model with single order parameter 
if and the potential Vq{}p) = Aip^ + Bip"^ + and we perturb it by adding 
Vi(0) = \eiip^diipdiip + h^e2ip^ + ■ We discuss the corresponding changes 
of the velocity of the planar interface. 
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1 Introduction 



Properties and evolution of interfaces is a very interesting topic in soft con- 
densed matter physics, see, e. g., [0 for a recent review. There are two main 
Unes of theoretical research in that direction. The first one, computer sim- 
ulations, provides a relation to microscopic, molecular dynamics level If. 
In the second approach one is satisfied with a coarse-grained description by 
effective Ginzburg-Landau models. 

In the framework of Ginzburg-Landau models the interfaces are described 
by particular solutions of pertinent nonlinear evolution equations. The non- 
linearity is crucial for existence of these solutions, but it is also a serious 
obstacle in obtaining them in an explicit form. In most cases only numerical 
solutions are available. A well-known exception is a model which appears in 
the context of react ion- diffusion chemical processes [^, and also is related to 
Landau-deGennes effective model for nematic liquid crystals f^. It involves 
a single, non-conserved, real- valued order parameter ip{x"',t), and the free 
energy functional has the form 

Fo = l d'x(^^\/^V^ + Voi^)y (1) 

Vo{ip)=A^^ + Bip'' + ip\ (2) 

where V = (gf^), ct = 1,2,3, are Cartesian coordinates in the space. 
For simplicity we have rescaled and Fq so that they are dimensionless 
and the coefficient in front of the quartic term in Vo{ip) is equal to 1. The 
constants A, B are positive, and A is restricted to the interval < A < 
Then V{ip) has two minima and a maximum between them. In this model 
exact formulas for profile and velocity of a planar interface are known, see, 
e.g., i. 

Exactly soluble models provide natural starting points for perturbative 
expansions. If the expression (1) for the free energy is changed a little bit, 
one may hope that the corresponding change of the planar interface can 
be calculated with the help of a perturbative expansion around the known 
solution. In this way we would obtain analytical results for velocity and 
profile of planar interfaces in a class of Ginzburg-Landau models. Moreover, 
this would open the way for investigations of evolution of curved interfaces 
in that class of models, because one could use the expansion in curvatures 
for which the knowledge of planar interfaces is the only prerequisite. 
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In the present paper we explore the idea of the perturbative expansion 
for the planar interfaces. We show that indeed such perturbative scheme can 
be constructed. Corrections to the profile are obtained in it as solutions of 
ordinary, linear differential equations, while corrections to the velocity follow 
from integrability conditions which are due to existence of zero-modes. This 
type of integrability conditions is well-known in statistical physics . As an 
example we consider perturbations of the form 



Values of the perturbation parameters 61,62,63 are such that Fq + 6F is 
bounded from below. The 61 correction to the gradient term introduces 
dependence on (p of the elastic constant K = 1 + 6iip'^. The 62 and 63 terms 
are corrections to the potential Vo{(p). The perturbations (3) are vanishing 
in the disordered phase [cp = 0), while in the ordered phase 7^ 0) they 
may give finite contributions. 

Our paper is organised as follows. In Section 2 we construct the pertur- 
bative expansion. In Section 3 we apply it to the model (1), (2) with the 
particular perturbation 6F given by formula (3). Summary and remarks are 
presented in Section 4. In the Appendix we show that the perturbatively 
calculated velocity of the planar interface coincides with the one calculated 
from an exact formula, provided that the perturbation is not too strong and 
that the perturbative series is convergent. 

2 General structure of the perturbative ex- 
pansion 

We will consider Ginzburg-Landau models in which time evolution of the 
order parameter ip is governed by equation of dissipative type, namely 



6F 




where 



Vi = -6iip^VipVip + -62^^ + ^£^3V5^- 



(3) 



dip 



(4) 



dt 
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where V = Vq + eVi, t denotes a rescaled dimensionless time, A = dxccd^oc is 
the three-dimensional Laplacian, and 



We have assumed that the perturbation Vi does not depend on second and 
higher derivatives of (/?. The Lh.s. of Eq. (4) is equal to the variational 
derivative dF/5ip multiplied by —1. 

The interface is identified with a solution of the evolution equation which 
interpolates between minima </?+, (p- of the potential V. Potential Vq given 
by formula (2) has two minima, and the perturbation Vi is by assumption 
weak enough, so that also V has just two minima which smoothly merge 
with the minima of Vq when Vi is switched off. In the present paper we 
concentrate on planar interfaces. Due to translational invariance of Eq.(4), 
it is sufficient to consider the planar interfaces parallel to the plane {x^,x'^) 
- then (/? is a function of and t. We require that in the limits x^ — > ±oo, 

(p{x^,t) (p±, d^sip 0. (5) 

Let us multiply Eq.(4), in which now V</7 —>■ d^i^^x^ ,t), by d^^^p and 
integrate it over x^ in the interval (— oo, +00). With the help of integration 
by parts, and using a formula for the full derivative of V with respect to x^, 
namely 

dV _ dVg ^ dV ^2 
dx^ d(f ^ ^ d{dx3(p) 
we obtain the following relation 

/■+00 
dx'dt^d,3ip = V{ip_)-Vi^+). (6) 
-00 

It implies that (p has to be time-dependent if the minima (p± are nondegen- 
erate. The simplest way to ensure that the l.h.s. of formula (6) is constant 
in time and nonvanishing consists in assuming that (p depends only on the 
combination x^ — vt, 

<p^<p{x^ -vt), (7) 

where v is a constant equal to the velocity of the interface. This form of 
dependence on time is also consistent with boundary conditions (5). For 
such (p, relation (6) can be written in the form 

np^)-nfA 

ii-dz{d,^r ' ^ ' 
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where z — — vt. Equation (4) is reduced to 

dl^ + vd,^-V'{^,d,^)=0, (9) 
where </? is a function of z only, and 

^{Z) ^ 99± (10) 

when z —>■ ±oo, respectively. The asymptotic values (fi± are determined from 
the equations 

V'{ip±) = 0. (11) 

From a mathematical viewpoint (8) and (9) should be regarded as a set of 
equations for v and (p{z), with (10) as the boundary conditions for (p{z). 
The perturbative Ansatz for ip{z) has the form 

ifi{z) = (po{z) + eipi{z) + e^(p2{z) + ■■■. (12) 

Here (poi^) is the initial, unperturbed interface. It obeys the following equa- 
tion 

d^ipo + vod^ipo - Vo{ipo) = 0, (13) 
with the boundary conditions at z — > ±oo 

cpo{z) a±, 

where a± are the minima of Vq. The velocity Vq of the unperturbed interface 
is given by the formula analogous to (8), namely 



Vo{a+) - Voja. ) 
j^^dz{d.^,y 



Formula (8) suggests that the velocity of the interface has the following per- 
turbative expansion 

V ^ Vq -\- £Vi -\- e^V2 H . (15) 

Let us stress that in equations (8), (9) and in expansion (12) the inde- 
pendent variable is just z, and not and t. Therefore, when solving Eqs. 
(8), (9), the variable z is not treated as — vt, and therefore we do not 
expand z in e, contrary to what formula (15) might suggest. Only after ip{z) 
and V are determined we may substitute z = x^ — vt in order to relate the 
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solution of the set of equations (8), (9) with the interface solution of the 
original evolution equation (4). 

Inserting the expansions (12), (15) in Eq.(9) and equating to zero the 
l.h.s. of (9) order by order in e, we obtain the infinite chain of equations for 
the corrections ^Pkik > 1. These equations have the form 

L(Pk = fk: (16) 

where 

L^dl + vod,-V^'{^o): (17) 

and fk in general depend on vi, . . . ,Vk and (fi,..., (fik-i- Simple calculations 
give 

/i = -vi d,(po + Vl{ipo, d,(po), (18) 

and for A; > 2 

fe-i 

fk = -Vkd^ifo - ^ Vj d^cpk-j + hk{z), (19) 

where hk{z) denotes the sum of all terms of the k-th order in e obtained 
by inserting formula (12) in Vq{(p) + eV({ip, dz^p) with exception of the term 
VQ{(po)'Pk which has been included into the l.h.s. of Eq.(16). 

It is a well-known fact that inhomogeneous linear equations not always 
have solutions - certain integrability conditions have to be satisfied. In the 
case of Eq.(16) such conditions appear because for the operator L there exists 
so called left zero- mode, that is a normalizable function ipi{z) such that the 
following identity holds 

/ + 00 
dziPi{z)Lipk = 0. (20) 
-oo 



If (fik obeys Eq.(16) then 



/+00 
-oo 



dz^i{z)fk = 0. (21) 



In our perturbative scheme the integrability conditions (21) serve as equa- 
tions which determine the perturbative contributions Vk to the velocity of 
the interface. Using formulas (18), (19) and conditions (21) we find that 

p+oo 

vi = / dzVliipo, d,ipo)M^), (22) 
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and for A; > 2 

fc-i 



Vk^N ^ j dzyhk{z)-Yj^ Vj d^ipk-j j M^)^ (23) 

where 



/+00 
dzi;r{z)M^)- 
-oo 



The functions ipijipr will be given shortly. In the Appendix we check that 
these recursive formulas and the expansion (15) give the velocity v which 
coincides with the one obtained by expanding in e the r.h.s. of formula (8), 
provided that e is not too large. 

The left zero-mode ipi can be found in the following way. Differentiation 
of the both sides of Eq.(13) with respect to z gives the identity 

L^r = 0, (24) 

where 

ijjr ^ d;,ipo(z). (25) 
The function t/jr is called the right zero-mode, t/ji has the form 

'ijji{z) ^ exjp{voz)'il}r{z). (26) 

Because derivation of identity (20) involves integration by parts, we have to 
discuss behaviour of ipi and ipr for z — > ±oo. Notice that Eq.(24) implies 
that for 2; ^ ±oo 



M^) = exp (-^^) exp + K"(a±) zj , 

respectively. Therefore ipi{z) exponentially vanishes for z ±00, provided 
that 

V^ia^) > 0. 

For Vq given by formula (2) this condition is satisfied if < A < as it 
has been assumed. When checking identity (20), boundary terms like, e.g., 
(dz'ilJi)ipk\l^, vanish because 

(fik ^ const for z —>■ ±00, (27) 
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in accordance with the boundary condition (10). 

Equations (16) can be solved in a standard way ||]. Adopting formulas 
given in 0] to the case at hand we find that 

= M^) r dCe^^^^MOfkiC) - Aiz) r rfCe"»^^2(C)/fc(C), (28) 

J-oo Jo 

where 

Mz)=M^) / rfCe-""^^,-2(C)- (29) 
Jo 

The functions ipr , V'2 form a pair of linearly independent solutions of the 
homogeneous equation Lip = 0. Notice that the integral present in the 
first term on the r.h.s. of formula (28) vanishes for z —>■ +oo due to the 
integrability conditions (21). It turns out that in the limits z —>■ ±oo both 
terms on the r.h.s. of formula (28) approach constants, and ifk obey the 
conditions (27). 

General solution of Eqs.(16) is obtained by adding to the r.h.s. of formula 
(28) the general solution of the corresponding homogeneous equation, that is 
the function Ckipr{,z) -\- dkip2{,z) . However, because 'ip2{z) exponentially grows 
for large |^| we have to put dk = 0. In order to determine the constants Ck 
we have to impose a condition on the interface solution ip{z) in addition to 
the boundary conditions (10). Actually, Eq.(9) and the boundary conditions 
(10) do not determine ip{z) uniquely - due to invariance with respect to 
translations in z we can take ip{z — zq) with arbitrary zq. We will require 
that 

<^(0)=<^o(0), (30) 

where (po{z) is a concrete, explicitly given function describing the unper- 
turbed interface with fixed location on the z axis. The condition (30) breaks 
the translational invariance. It implies that 

MO) = 0. (31) 

The solutions (pk{z) given by formula (28) obey this condition, therefore also 
the constants Ck are vanishing. Thus, the perturbative scheme supplemented 
with the condition (30) yields unique interface solution (p{z). In this solution 
we can of course shift the variable z, that is to substitute z ^ z — Zq simul- 
taneously in all contributions (po,(pk, in order to obtain the other solutions 
implied by the translational invariance. 
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3 Effects of the perturbations of the form 

Vi = \eiLp'^V LpV LP + \e2(f^ + ^£3^. 

Let us start from a description of the unperturbed interface ^o{z). The 
evolution equation has the form (4) with V replaced by Vq given by formula 
(2). The interface solution is well-known, see, e.g., [Q. Let us quote the 
relevant formulas from and [^]. We use the following abbreviations 

z 1 



4/o' 2v^/o' 
/o^^ = ^ (35 + V9fi2 _ 32A) . (32) 



where 

Then 

and 



V5o(^) = ^ (1 + tanhs) , (33) 



a+ 1 



oLq cosh S 
2/ ^1 /""^ 

^2(z) = ^ / rfaexp[-2(2a + l)a] [1 + exp(2(T) 

fl-i- cosh s Jo 



cosh s Jo 

The integral over a can be easily calculated, but the resulting formula for 1^2 
is quite long. In the absence of perturbations, z = — vot, where 

q 

vo = ^ [b - _ 32^^ . (34) 

The interface exists and it is stable when < A < Usually, the 

coefficient B is regarded as independent of temperature, while A = a(T — T^), 
where a > 0. The inequality given above fixes the temperature range (T^,, Tc) 
in which the interface exists. Notice that at the temperature Tq such that 
A = ^ the velocity vq vanishes. It is clear that T* < Tq < Tc. The left 
zero-mode has the form 

^Pliz) = ^exp[2(2a- 1).]^^, (35) 
o<o cosh s 
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where ^ 

The value of the parameter a monotonically increases from to 1 when A 
varies from to 

The first order correction v\ to the velocity is calculated from formula 
(22). In the present case 

^(V^o) = -eiV'oV'r^ - eiV^o^^^^r + £2V5o^ + ^sV^o^- (36) 

The integrations over z ( or equivalently over s ) can be related to the Euler 
B function 

/•oo t^~^ 

B(x,y) = / dt- ^— , 

by the substitution y = exps. The resulting expressions can be simplified 
with the help of the well-known formulas 

r(x)r(v) 

B{x,y)= ) , T{x + l)=xT{x). 

T{x + y) 

After straightforward calculations we obtain 



Vl 



(a + l)(2a + l) 



160/, 
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L /o 2a + 3 (2a + 3)(a + 2) 

[a - -)6i - 62 + — —7- ^ €3 

2 y/2 1 — a 12(1 — a) 



(37) 

Formula (37) shows that the 62, £3 corrections become more and more 
pronounced as a increases towards 1, that is as the temperature increases 
towards Tc. The singularity at a = 1 is due to the fact that for a = 1 the 
left zero-mode tpi and VY(v^o) do not vanish when s — > +00, and therefore 
the integral in formula (22) is divergent. 

The El term in formula (37) vanishes for a = \. The reason is that the 
corresponding term in the perturbation potential Vi does not contribute to 
the values of the full potential V at the two minima ip±. Therefore, it can 
influence the velocity only by changing the denominator in formula (8). For 
a = I the numerator in (8) vanishes if £2 = ^3 = and the £1 term can not 
contribute to the velocity. 

The facts that 62 term in (37) is negative, and that the £3 term is posi- 
tive, can be explained by the observation that the 62 term in the potential 
diminishes the potential energy difference across the interface, while the 
term in the potential works in the opposite direction. 
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4 Summary and discussion 

1. The purpose of this work is to provide a practical tool for analytic, ap- 
proximate computations of characteristics of the perturbed planar interfaces 
within the framework of Ginzburg-Landau models. We have shown how one 
can systematically calculate corrections to the profile and the velocity of 
the interface. In particular, the scheme can be applied to perturbations of 
domain walls, which can be regarded as special, static interfaces appearing 
when the two minima of Vq are degenerate. In the example discussed in Sec- 
tion 3 this is the case when a = 1/2. The scheme is relatively simple. The 
final formulas (22), (23), (28), (29) contain one dimensional integrals which 
always can be tackled numerically, and in many cases calculated analytically. 

2. We have not delved into the problem of convergence of the perturbative 
series. Such a mathematical investigation does not belong to the scope of 
our work. It is clear that if the perturbations shift the two minima of Vq only 
slightly, the corresponding changes of the profile and of the velocity of the in- 
terface are also small, hence they can be calculated as small corrections. This 
does not necessarily mean that the series is convergent - it can belong to a 
wider class of asymptotic series, but it is acceptable from physical viewpoint. 

3. Using the results of the present paper in combination with the curva- 
ture expansion mentioned in the Introduction, one can calculate evolution of 
curved interfaces in the perturbed Ginzburg-Landau models. Comparisons 
of such theoretical predictions with experimental observations of evolution 
of the interfaces could help to determine the best formula for the free en- 
ergy of the system. Let us remind that the formulas for the free energy in 
the Ginzburg-Landau models rarely can be derived from underlying micro- 
scopic theories. In most cases they are postulated on basis of qualitative 
phenomenological considerations. General discussion of various formulas for 
the free energy in the context of theory of liquid crystals can be found in 

Hi- 

5 Appendix. Resummation of the perturba- 
tive contributions to the velocity 

The perturbative corrections Vk to the velocity of the interface are given by 
formulas (22), (23). They have been obtained in somewhat indirect way, 
namely from the integrability conditions. On the other hand, the velocity is 
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given by the exact formula (8). We would like to check that the sum (15) of 
the perturbative contributions Vk, below denoted by v, 

oo 

V ^vo + + £^^k. (38) 



k=2 



does coincide with the velocity v given by formula (8). Actually, because for- 
mula (8) follows directly from Eq.(9), it is sufficient to check that v coincides 
with the V present in that equation. 

In the first step, wc substitute for Vq, Vi and Vk { k > 2 ) in (38) formulas 
(14), (22) and (23), respectively. Next, we notice that the definition of hk 
given below formula (19) is equivalent to the following formula 

oo 
k=2 

where 

Vo{(p) + eVl{(p, d^if) = d^^if + vd^if, 

and 

according to Eqs. (9), (13) and definition (17). We also use the formula 

oo k—1 / oo \ / oo \ 

k=2 j=l \i=l / \j=l J 

After a simple calculation we obtain from formula (38) that 

/+00 
dzi/jid^ifi = 0. (39) 
-oo 



V = V, 



Thus, indeed 
provided that 

Because 

dz(p = V'r + edz(pi H , 

that last condition is satisfied if e is small enough, that is if the perturbation 
is not too strong. 



/+0O 
-oo 
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